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Abstract 

A bound state in the continuum (BIC) is an unusual localized state that is embedded in a 
continuum of extended states. Here, we present the general condition for BICs to arise from 
wave equation separability and show that the directionality and dimensionality of their resonant 
radiation can be controlled by exploiting perturbations of certain symmetry. Using this general 
framework, we construct new examples of separable BICs in realistic models of optical potentials 
for ultracold atoms, photonic systems, and systems described by tight binding. Such BICs with 
easily reconfigurable radiation patterns allow for applications such as the storage and release of 
waves at a controllable rate and direction, systems that switch between different dimensions of 
confinement, and experimental realizations in atomic, optical, and electronic systems. 
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I. INTRODUCTION 


The usual frequency spectrum of waves in an inhomogeneous medium consists of localized 
waves whose frequencies he outside a continuum of delocalized propagating waves. For some 
rare systems, the bound states can be at the same frequency as an extended state. Such 
bound states in the continuum (BICs) represent a method of wave localization that does 
not forbid outgoing waves from propagating, unlike traditional methods of localizing waves 
such as potential wells in quantum mechanics, conducting mirrors in optics, band-gaps, and 
Anderson localization. 

BICs were first predicted theoretically in quantum mechanics by von Neumann and 
Wigner [1]. However, their BIC-supporting potential was highly oscillatory and could not 
be implemented in reality. The BIC was seen as a mathematical curiosity for the following 
50 years until Friedrich and Wintgen showed that BICs could exist in models of multi¬ 
electron atoms |2]. Since then, many theoretical examples of BICs have been demonstrated 
in quantum mechanics M, electromagnetism [TnH23], acoustics (TU], and water waves 
with some experimentally realized [251 - I5U] . A few mechanisms explain most examples of 
BICs that have been discovered. The simplest mechanism is symmetry protection, in 
which the BIC is decoupled from its degenerate continuua by symmetry mismatch mu 
[ini [ig [251 |26l [31]. Or, when two resonances are coupled, they can form BICs through 
Fabry-Perot interference [SI [IHl [23 122] or through radiative coupling with tuned resonance 
frequencies [21 [231 ED]. More generally, parametric BICs can occur for special values of 
parameters in the Hamiltonian [Ml HZl [2T| [22| [271132] , which can in some cases be understood 
using topological arguments [33] . 

There exists yet another class of BICs, which has thus far been barely explored, and 
literature covering it is rare and scattered; these are BICs due to separability [3H6l[nHI3|. 
In this paper, we develop the most general properties of separable BICs by providing the 
general criteria for their existence and by characterizing the continuua of these BICs. We 
demonstrate that separable BICs enable control over both the directionality and the dimen¬ 
sionality of resonantly emitted waves by exploiting symmetry, which is not possible in other 
classes of BICs. Our findings may lead to applications such as tunable-Q directional res¬ 
onators and quantum systems which can be switched between quantum dot (OD), quantum 
wire (ID), and quantum well (2D) modes of operation. While previous works on separable 
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BICs have been confined to somewhat artificial systems, we also present readily experimen¬ 
tally observable examples of separable BICs in photonic systems with directionally resonant 
states, as well as examples in optical potentials for cold atoms. 


II. RESULTS 

A. General Condition for Separable BICs 

We start by considering a simple example of a two-dimensional separable system - one 
where the Hamiltonian operator can be written as a sum of Hamiltonians that each act on 
distinct variables x or y, i.e; 

H = H,{x) + Hy{y). (1) 

Denoting the eigenstates of and Hy as 'ipl.{x) and ip^iy), with energies E^. and 
it follows that H is diagonalized by the basis of product states with energies 

= El + E^. If Hx and Hy each have a continuum of extended states starting at zero 
energy, and these Hamiltonians each have at least one bound state, and respectively, 
then the continuum of H starts at energy min(i?°, Ey), where the 0 subscript denotes ground 
states. Therefore, if there exists a bound state satisfying E^’^ > min(i?°, i?°), then it is a 
bound state in the continuum. We schematically illustrate this condition being satisfied 
in Figure 1, where we illustrate a separable system in which Hx has two bound states, 
and ipl, and Hy has one bound state, ipy. The first excited state of Hx combined with the 
ground state of Hy, iplipy, while spatially bounded, has larger energy, El + Ey, than the 
lowest continuum energy, E^., and is therefore a BIC. 

To extend separability to a Hamiltonian with a larger number of separable degrees of free¬ 
dom is then straightforward. The Hamiltonian can be expressed in tensor product notation 

N 

as if = ^ Hi, where Hi = ®hi® In this expression, N is the number of sepa- 

i=l 

rated degrees of freedom, hi is the operator acting on the Tth variable, and I is the identity 
operator. The variables may refer to the different particles in a non-interacting multi-particle 
system, or the spatial and polarization degrees of freedom of a single-particle system. De¬ 
note the rijih eigenstate of hi by with energy i?”^. Then, the overall Hamiltonian H is 
trivially diagonalized by the product states |ni, n 2 , • • •, with 

corresponding energies E = E^^ + E^^ 4- ■ ■ ■ Epp' . Denoting the ground state of hi by E^ 
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and defining the zeros of the hi such that their continuua of extended states start at energy 
zero, the continuum of the overall Hamiltonian starts at min iE^). Then, if the separated 

l<i<N * 

operators {hi} are such that there exists a combination of separated bound states satisfying 

N 

E = Ey > min (Ef), this combined bound state is a BIC of the overall Hamiltonian. 

i=l l<i<Af ^ 

For such separable BICs, coupling to the continuum is forbidden by the separability of the 
Hamiltonian. Stating the most general criteria for separable BICs allows us to straight¬ 
forwardly extend the known separable BICs to systems in three dimensions, multi-particle 
systems, and systems described by the tight-binding approximation. 

B. Properties of the Degenerate Continuua 

A unique property that holds for all separable BICs is that the delocalized modes de¬ 
generate to the BIC are always guided in at least one direction. In many cases, there are 
multiple degenerate delocalized modes and they are guided in different directions. In 2D, 
we can associate partial widths Fa, and F^, to the resulting resonances of these systems when 
the separable system is perturbed. These partial widths are the decay rates associated with 
energy-conserving transitions between the BIC and states delocalized in the x and y direc¬ 
tions under perturbation, respectively. When separability is broken, we generally can not 
decouple leakage in the x and y directions because the purely-a; and purely-y delocalized 
continuum states mix. However, in this section, we show that one can control the radiation 
to be towards the x (or y) directions only, by exploiting the symmetry of the perturbation. 

In Fig. 2(a), we show a two-dimensional potential in a Schrodinger-like equation which 
supports a separable BIC. This potential is a sum of Gaussian wells in the x and y-directions, 
given by 

_2cc 2 _2ld 

V{x,y) = -I4e ^ - VyC . (2) 

This type of potential is representative of realistic optical potentials for ultracold atoms 
and potentials in photonic systems. Solving the time-independent Schrodinger equation for 
{14, By} = {1.4, 2.2} and {axiCXy} = {5,4} (in arbitrary units) gives the spectra shown in 
Fig. 2(b). Due to the x and y mirror symmetries of the system, the modes have either even 
or odd parity in both the x and the y directions. This system has several BICs. Here we 
focus on the BIC \nx-,ny) = |2,1) at energy E"^'^ = —1.04, with the mode prohle shown in 
Fig. 2(c); being the second excited state in x and the hrst excited state in y, this BIC is 


4 



even in x and odd in y. It is only degenerate to continuum modes |0, — E^) extended in 

the y direction (Fig. 2(d)), and \E‘^’^ — E^, 0) extended in the x direction (Fig. 2(e)), where 
an E label inside a ket denotes an extended state with energy E. 

If we choose a perturbation that preserves the mirror symmetry in the y-direction, as 
shown in the inset of Fig. 2(f), then the perturbed system still exhibits mirror symmetry in 
y. Since the BIC |2,1) is odd in y and yet the x-delocalized continuum states \E‘^'^ — i?°,0) 
are even in y, there is no coupling between the two. As a result, the perturbed state radiates 
only in the y direction, as shown in the calculated mode prohle in Fig. 2(f). This directional 
coupling is a result of symmetry, and so it holds for arbitrary perturbation strengths. 

On the other hand, if we apply a perturbation that is odd in x but not even in y, as shown 
in the inset of Fig. 2(g), then there is radiation in the x direction only to hrst-order in time- 
dependent perturbation theory. Specihcally, for weak perturbations of the Hamiltonian, 5V, 
the hrst-order leakage rate is given by Fermi’s Golden Rule for bound-to-continuum coupling, 
F ~ 1(2, 1|5 R|'^c)PPc(A'^’^); where Pc{E) is the density of states of continuum c, and c 

labels the distinct continuua which have states at the same energy as the BIC. Since the 
BIC and the y-delocalized continuum states |0, E"^’^ — E^) are both even in x, the odd-in-x 
perturbation does not couple the two modes directly, and Fy is zero to the hrst order. As 
a result, the perturbed state radiates only in the x direction, as shown in Fig. 2(g). At 
the second order in time-dependent perturbation theory, the BIC can make transitions to 
intermediate states k at any energy, and thus the second-order transition rate, proportional 
to = Xlc I Xlfe vanish because the intermediate state can 

have even parity in the x-direction. 

Another unique aspect of separability is that by using separable BICs in 3D, the di¬ 
mensionality of the conhnement of a wave can be switched between one, two, and three by 
tailoring perturbations applied to a single BIC mode. The ability to do this allows for a 
device which can simultaneously act as a quantum well, a quantum wire, and a quantum 
dot. We demonstrate this degree of control using a separable potential generated by the 
sum of three Gaussian wells (in x, y, and directions) of the form in Equation 2 with 
strengths {0.4, 0.4,1} and widths {12,12,3}, all in arbitrary units. The identical x and y 
potentials have four bound states at energies E^ = Ey = —0.33, —0.20, —0.10 and — 0.029. 
The potential has two bound states at energies E^ = —0.61 and — 0.059. All of the 
continuua start at zero. The BIC state |1,1,1) at energy = —0.47 is degenerate to 
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|0, 0, Ez), |0,1, Ez), |1, 0, Ez), \Ex, m, 0), \n, Ey, 0) and \Ez;, Ey, 0), where Ei denotes an energy- 
above zero, and m and n denote bonnd states of the x and y wells. For pertnrbations which 
are even in the z direction, the BIG |1,1,1) does not conple to states delocalized in x and 
y {\Ex, m, 0), \n, Ey, 0), and lE^, Ey, 0)), becanse they have opposite parity in z, so the BIG 
radiates in the 2 ;-direction only. On the other hand, for pertnrbations which are even in the 
X and y directions, the conpling to states delocalized in z (|0,0, Ez), |0,1, Ez), and |1, 0, Ez)) 
vanishes, meaning that the BIG radiates in the xy-plane. 


C. Proposals for the experimental realizations of separable BICs 


BIGs are generally difficnlt to experimentally realize becanse they are fragile nnder per¬ 
tnrbations of system parameters. On the other hand, separable BIGs are straightforward to 
realize and also robust with respect to changes in parameters that preserve the separability 
of the system. This reduces the difficulty of experimentally realizing separable BIGs to en¬ 
suring separability. And while this is still in general nontrivial, it can be straightforwardly 
achieved in systems that respond linearly to the intensity of electromagnetic helds. In the 
next two realistic examples of BIGs that we demonstrate, we use detuned light sheets to 
generate separable potentials in photonic systems and also for ultracold atoms. 

Paraxial optical systems — As a hrst example of this, consider electromagnetic waves prop¬ 
agating paraxially along the 2 ;-direction, in an optical medium with spatially non-uniform 
index of refraction n{x,y) =710-1- ^n{x,y), where rio is the constant background index of 
refraction and 6n -C hq. The slowly varying amplitude of the electric held '0(a;,i/, z) satishes 
the two-dimensional Schrodinger equation (see Ref. 34 and references therein). 




k6n , 

—V’- 

no 


(3) 


Here, = d‘^/dx‘^ + d‘^/dy‘^, and k = 27rno/A, where A is the wavelength in vacuum. 
The modes of the potential 5n[x, y) are of the form ip = Aj{x, where Aj is the prohle, 

and Pj the propagation constant of the jth mode. In the simulations we use no = 2.3, and 
A = 485 nm. 

Experimentally, there are several ways of producing potentials 6n{x, y) of different type, 
e.g., periodic [54], random [35], and quasicrystal [53]. One of the very useful techniques is the 
so-called optical induction technique where the potential is generated in a photosensitive 
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material (e.g., photorefractives) by employing light [S^. We consider here a potential gen¬ 
erated by two perpendicular light sheets, which are slightly detuned in frequency, such that 
the time-averaged interference vanishes and the total intensity is the sum of the intensities 
of the individual light sheets. The light sheets are much narrower in one dimension {x for 
one, y for the other) than the other two, and therefore each light sheet can be approximated 
as having an intensity that depends only on one coordinate, making the index contrast 
separable. This is schematically illustrated in Fig. 3(a). 

If the sheets are Gaussian along the narrow dimension, then the potential is of the form 
6 n{x,y) = 5no[exp(—2(x/(j)^) -|-exp(—2(?//cr)^)]. It is reasonable to use a = 30 pm and 
6no = 5.7 X 10“^. For these parameters, the one dimensional Gaussian wells have four 
bound states, with f3 values of (in mm“^): 2.1,1.3,0.55, and 0.11. There are eight BIGs: 
|1, 2), |2,1), |1, 3), |3,1), |2,3), |3,2),|2,2), and |3,3). Among them, |1,3) and |3,1) are 
symmetry protected. Additionally, the BIGs |1,3) and |3,1) can be used to demonstrate 
directional resonance in the x oi y directions, respectively, by applying perturbations even 
in X or p, respectively. Therefore, this photonic system serves as a platform to demonstrate 
both separable BIGs and directional resonances. 

Optical potentials for ultracold atoms — The next example that we consider can serve as 
a platform for the hrst experimental measurement of BIGs in quantum mechanics. Gonsider 
a non-interacting neutral Bose gas in an optical potential. Optical potentials are created by 
employing light sufficiently detuned from the resonance frequencies of the atom, where the 
scattering due to spontaneous emission can be neglected, and the atoms are approximated 
as moving in a conservative potential. The macroscopic wavefunction of the system is then 
determined by solving the Schrodinger equation with a potential that is proportional to the 
intensity of the light [38] . 

As an explicit example, consider an ultracold Bose gas of atoms. An optical potential 
is made by three Gaussian light sheets with equal intensity R = R = I 3 , widths 20 pm, 
and wavelengths centered at A = 1064 nm [39]. This is schematically illustrated in Fig. 
3(b). The intensity is such that this potential has depth equal to ten times the recoil energy 
Er = 2 mx 2 ■ By solving the Schrodinger equation numerically, we hnd many BIGs; the 
continuum energy starts at a reduced energy R = = —296.24, where Xq is chosen 

to be 1 fim. The reduced depth of the trap is —446.93. Each one-dimensional Gaussian 
supports 138 bound states. There are very many BIGs in such a system. For concreteness. 
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an example of one is |30,96,96), with reduced energy —146.62. 

Tight Binding Models— The hnal example that we consider here is an extension of the 
separable BIC formalism to systems which are well-approximated by a tight-binding Hamil¬ 
tonian that is separable. Consider the following one-dimensional tight-binding Hamiltonian, 
Hi, which models a one-dimensional lattice of non-identical sites: 

Hi = '^e\\k){k\+ti^{\l){m\ + |m)(/|), (4) 

k (Im) 

where {Im) denotes nearest-neighbors, ef is the on-site energy of site k, and k, I, and m run 
from —cx) to ex. Suppose ef = —V for \k\ < N, and zero otherwise. For two Hamiltonians 
of this form, Hi and H 2 , Hi ® I + I ® H 2 describes the lattice in Fig. 3(c). If we take 
Hi = H 2 with {V,t,N} = {—1, —0.3,2}, in arbitrary units, the bound state energies of the 
ID-lattices are numerically determined to be —0.93,-0.74,-0.46, and — 0.16. Therefore 
the states |2, 2), |2, 3), |3,2), |3, 3), |3,1) and |1,3) are BICs. The last two of these are also 
symmetry-protected from the continuum as they are odd in x and y while the four degenerate 
continuum states are always even in at least one direction. Of course, many different physical 
systems can be adjusted to approximate the system from Eq. (4), so this opens a path for 
observing separable BICs in a wide variety of systems. 

III. SUMMARY 

With the general criterion for separable BICs, we have extended the existing handful of 
examples to a wide variety of wave systems including: three-dimensional quantum mechan¬ 
ics, paraxial optics, and lattice models which can describe 2D waveguide arrays, quantum dot 
arrays, optical lattices, and solids. These BICs exist in other wave equations such as the 2D 
Maxwell equations, and the inhomogeneous scalar and vector Helmholtz equations, meaning 
that their physical domain of existence extends to 2D photonic crystals, microwave optics, 
and also acoustics. To our knowledge, they do not exist (except perhaps approximately) in 
isotropic media in the full 3D Maxwell equations VxVxE = V(V-E) — V^E = e^E 
because the operator V(V ■ ) couples the different coordinates and renders the equation 
non-separable. We have demonstrated numerically the existence of separable BICs in mod¬ 
els of trapped Bose-Einstein condensates and photonic potentials, showing that separable 
BICs can be found in realistic systems. These simple and realistic models can facilitate the 



observation of BICs in quantum mechanics, which to this date has not been conclusively 
done. 

More importantly, we have demonstrated two new properties unique to separable BICs: 
the ability to control the direction of emission of the resonance using perturbations, and 
also the ability to control the dimensionality of the resulting resonance. This may lead 
to two applications. In the hrst, perturbations are used as a switch which can couple 
waves into a cavity, store them, and release them in a hxed direction. In the second, the 
number of dimensions of conhnement of a wave can be switched between one, two, and three 
by exploiting perturbation parity. The property of dimensional and directional control of 
resonant radiation serves as another potential advantage of BICs over traditional methods 
of localization. 
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Hx(x) + Hy(y) = H(x,y) 
E 



FIG. 1: A schematic illustration demonstrating the concept of a separable BIG in two dimensions. 
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FIG. 2: (a) A separable potential which is a sum of a purely x-dependent Gaussian well and a 
purely y-dependent Gaussian well, (b) The relevant states of the spectrum of the x-potential,?/- 
potential, and total potential, (c) A BIC supported by this double well. (d,e) Continuum states 
degenerate in energy to the BIC. (f) A y-delocalized continuum state resulting from an even- 
y-parity perturbation of the BIC supporting potential, (g). An x-delocalized continuum state 
resulting from an odd-x-parity perturbation. 
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FIG. 3: Separable physical systems with BICs. (a) A photorefractive optical crystal whose index is 
weakly modified by two detuned intersecting light sheets with different intensities, (b) An optical 
potential formed by the intersection of three slightly detuned light sheets with different intensities, 
(c) A tight-binding lattice. 
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